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A mirror scattering vacuum fields is submitted to a quan- 
tum fluctuating radiation pressure. It also experiences a mo- 
tional force, related to force fluctuations through fluctuation- 
dissipation relations. The resulting position fluctuations of 
the coupled mirror are related to the dissipative part of the 
mechanical admittance. We compute the time dependent po- 
sition commutator, which makes apparent the difference be- 
tween the low-frequency and high-frequency masses, and the 
anticommutator noise spectrum, which describes the ultimate 
sensitivity in a length measurement using mirrors. 
PACS: 03.65 - 12.20 - 42.50 



INTRODUCTION 

A mirror which scatters a laser field is submitted to 
a fluctuating radiation pressure Q, which is associated 
through fluctuation-dissipation relations with a motional 
force, that is a mean force for a moving mirror Q|. The 
resulting random motion has been studied in detail since 
it determines the ultimate sensitivity of the interferome- 
ters designed for gravitational wave detection || . 

In the absence of laser irradiation, the mirror still scat- 
ters vacuum fields, so that the radiation pressure fluctu- 
ates Pj, and the motional force does not vanish [§,1). 
Fluctuations and dissipation are also directly connected 
in this case 0. As a consequence of its coupling with 
vacuum radiation pressure, the mirror's position has to 
acquire some quantum fluctuations, even if it is intro- 
duced as a classical variable. On the other hand, the 
mirror has proper quantum fluctuations, even if it is not 
coupled to the field. The purpose of the present paper 
is to elucidate the relation between those acquired and 
proper fluctuations. 

In order to reach consistent results, it is necessary that 
the mirror behave as a stable system when coupled to vac- 
uum. For a perfect mirror however, 'runaway solutions' 
appear, which are analogous to the ones encountered in 
classical electron theory || . We will use the fact that this 
instability problem is solved M by considering a partially 
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transmitting mirror, causal and transparent at frequen- 
cies higher than a reflection cutoff u>c much smaller in 
reduced units than the mirror's mass mo 

huJc -C m c 2 (1) 

This condition, which allows one to ignore the recoil ef- 
fect in the derivation of the motional force, also implies 
passivity, that is the incapacity of vacuum to sustain run- 
away solutions. Then, stability results from passivity [0. 

We shall obtain that the position fluctuations of the 
coupled mirror are connected to the dissipative part of 
the mechanical admittance, in analogy with the thermal- 
ized position of an atomic system scattering a thermal 
field jlO| , or the thermalized voltage of an electrical os- 
cillator coupled to a resistor at thermal equilibrium [ pr| . 
However, at zero temperature, this equilibrium involves 
vacuum fluctuations rather than thermal ones p^ j. The 
non-commutativity of the force correlation functions will 
play a central role in this paper, as well as the relation of 
the commutator with the dissipative part of the motional 
susceptibility Jig ]. 

First, we study the 'input fluctuations', i.e. the fluc- 
tuations when coupling is disregarded. The mirror's dy- 
namics are described by two noise spectra correspond- 
ing respectively to the commutator and anticommutator, 
which are related to a mechanical admittance function 
and which behave resonantly in the vicinity of suspension 
eigenfrequencies. When coupling to vacuum radiation 
pressure is taken into consideration, the admittance func- 
tion is modified. Using the techniques of linear response 
theory jl3|], we deduce coupled fluctuations which are 
related to the modified admittance through fluctuation- 
dissipation relations. In this derivation, input position 
fluctuations are considered on an equal foot with input 
force fluctuations. Yet, it appears that the fluctuations of 
the coupled variables are completely determined by those 
of the input force. This means that the results might be 
obtained from a Langevin equation Jl4| , where all fluctu- 
ations are fed by the input force fluctuations; this simpler 
derivation is presented in Appendix A. 

Then, we discuss the two noise spectra related respec- 
tively to the commutator and to the anticommutator for 
the coupled position. The former describes a modifica- 
tion of the canonical commutation relation. Its time de- 
pendence, analyzed in Appendix B, appears to be con- 
nected to a difference between the low-frequency and 
high-frequency values of mass. The latter gives the quan- 
tum noise on the position of the coupled mirror. When 
integrated over frequency, it provides the position vari- 
ance; as a function of time, it describes the quantum 
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diffusion of a mirror coupled to vacuum radiation pres- 
sure. The anticommutator noise spectrum is involved in 
ultimate quantum limits for a position measurement || . 

We also analyse the autocorrelation of the coupled 
force and its cross correlation with the coupled position. 

For the sake of simplicity, our attention will focus upon 
the simple problem of an harmonically bound mirror. 
However, linear response formalism can also be used for 
the anharmonic oscillator, as shown in Appendix C. 



INPUT POSITION FLUCTUATIONS 

In a stationary state, we will write any correlation 
function as the sum of an antisymmetric part, associated 
with a commutator, and of a symmetric one, related to 
an anticommutator 



One also defines a mechanical susceptibility which 
describes the linear response of the position to an exerted 
force Ol 
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XqqM = 



m (luq - (uj + ie) 2 ) 



(3g) 



The parameter e — > + is inserted in order to ensure the 
causal character of the susceptibility. The susceptibility 
is directly connected to the mechanical impedance Z ln or 
to the mechanical admittance Y m 



i 



(4) 



Those functions obey the fluctuation-dissipation rela- 
tions |l3j applied to a quantum system at zero tempera- 
ture 



Cju.it) = (A(t)A(0)) - (A) 2 = h (a AA (t) + UA(t)) (2a) 

C AA {t) - C A A(-t) = (A(t)A(0) - A(0)A(t)) 

= 2hU A (t) (2b) 

C AA {t) + C AA (-t) = (A(t)A(0) - A(0)A(t)) - 2 (A) 2 
= 2ha AA (t) (2c) 

The antisymmetric part £ AA , typical of quantum fluctua- 
tions, would vanish for any classical stochastic treatment. 
We will denote for any function / 



The function £aaM is an odd function of uj while the 
noise spectrum <r AA [uj] is an even one. If A is a hermitean 
operator, both are real functions of uj. 

Let us first express the position fluctuations of an 
harmonic oscillator (mass mo, eigenfrequency luq) in its 
ground state within the linear response formalism. The 
correlation functions of this 'input position' q ln are easily 
computed (we use shorthand notation where stands 

for tqin q m) 



qq ^ 2mQUJo 
^ qq ^ ^ 2ijioujq 



qq 2mQUJo 
as well as the associated spectra 
Tin 



exp(-iujot) 

sin(cJo^) 
cos(ujot) 



-8(ui - cj ) 



m uj 

4" H = ^j- q (S(l> - ujo) - S(u> + u; )) 



qq 2rriQUJo 



(5(u) - uj ) + S(uj + uj j) 



(3a) 
(3b) 
(3c) 

(3d) 
(3e) 
(3f) 



C q \[w] = 2M{uj)Q[u>] (5a) 
H = M e(w) = 0(w)- 9(-uj) (5b) 

2i££ M = X* [w] - X? q [-w] = 2 l Im X f q [w] (5c) 

The 9(uj) function may be considered as the limit at zero 
temperature of Planck spectrum |lq| 
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Quantities evaluated from the correlation functions at 
equal times, like the canonical commutator or the vari- 
ances, may be expressed as integrals over frequency of the 
spectra £gg[a/] an d fggM- For instance, the commutator 
of two positions evaluated at different times 

<[«(*),«(0)]> = 2^ w (t) 

is a real and odd function of t. The commutator between 
velocity and position is simply deduced (v(t) is the mir- 
ror's velocity q'(t)) 



{[v(t),q(0)})=2ht vq (t) = 2K„(t) 



(6) 



Introducing a specific notation for integrals over fre- 
quency 
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(111.' 
2^ 



ffM = 5(0) 



the equal-time canonical commutator is 

([«(0), 9 (0)]) = 2^,(0) =-2*n^ 



(7) 



(8) 



The symmetrical correlation function (we assume that 
the mean value (q) vanishes) 

(q(t)q(0) + q(0)q(t)) = 2Tia qq (t) 

is a real and even function of t, and the stationary vari- 
ances are 
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Aq 2 
Av 2 



<?(0) 2 ) 

(V(0) 2 ) = hU^g 



(9a) 
(9b) 



Using equations (3) , one recovers the well known expres- 
sions 



(lv(0), q (0)}) = — 
m 



Aq 2 



Av 2 



2mo 



(10) 
(11) 



The covariance between position and velocity vanishes 
because cr q9 [w] is an even function of u>. 

The foregoing description of quantum fluctuations may 
seem too sophisticated for such a simple system as an 
harmonic oscillator. We see in the following that this 
formalism is well suited to the analysis of fluctuations 
for the system coupled to vacuum. Relations (5,6,8,9) 
will still be valid, with different expressions for the noise 
spectra, and will have a significant content in this case. 

Furthermore, expressions (5) for the spectra contain 
more information than the integrated quantities, even for 
the uncoupled variables. Consider for instance that po- 
sition is measured and that the measurement frequency 
band is characterized by a function G[u>] with a maxi- 
mum value of 1 at the frequency of the expected signal. 
Assuming that quantum noise may be dealed with in the 
same manner as thermal noise, one gets the effective noise 
as an integral over the detection bandwidth (with the no- 
tation of eq. 7) 



HGa, 



<i'i 



This effective noise Aq^ may be smaller than the position 
variance Aq 2 . It even vanishes if the suspension eigen- 
frequencies are outside the detection bandwidth, which is 
the case in any high sensitivity measurement; in the op- 
posite case, the noise would be of the order of the position 
variance Aq 2 . This explains why the proper quantum 
fluctuations of the mirror's position may be neglected 
when analyzing the limits in position measurements. At 
the end of the paper, we will come back to this discussion 
and show that the dissipative coupling between mirror 
and vacuum sets an ultimate quantum limit in position 
measurements. 



INPUT FORCE FLUCTUATIONS 

We now consider the input force fluctuations F m , that 
is the fluctuations of vacuum radiation pressure com- 
puted when disregarding reaction of the mirror's motion. 

These fluctuations are characterized by stationary cor- 



relation functions denoted Cp F , £p F and <7™ F , which 
are connected through fluctuation-dissipation relations 
to the linear susceptibility Xff describing the motional 
force for small displacements 



2^M = AH-x^W = 



(12a) 
(12b) 

2i Imx'^H (12c) 

These relations, which can be regarded as consequences 
of linear response theory, have been directly checked for 
the problem of a mirror coupled to vacuum radiation 
pressure |7]j|]. In particular, equations (12a, b) corre- 
spond to the fact that vacuum is the equilibrium state at 
zero temperature, so that radiation pressure can damp 
the mirror's motion, but cannot excite it. 

For a perfect mirror scattering a scalar vacuum field in 
a two-dimensional (2D) spacetime, the motional force is 
proportional to the third derivative of the position and 
corresponds to the following susceptibility, obtained by 
linearizing the expressions obtained for a perfect mirror 
in reference ||, or by considering the limit of perfect 
reflection in reference j?]] 



XffM 



ihuj 3 
6ttc 2 



The time constant r characterizes the weak coupling of 
mirror to vacuum 

h 



Girmoc 2 



A partially transmitting mirror ]7j is described by 
causal reflectivity and transmittivity functions with a re- 
flection cutoff. This approach finds many justifications. 
It leads to an expression for the mean Casimir force be- 
tween two mirrors, which is free from the divergences 
usually associated with the infiniteness of vacuum energy 
|]l6f . It also provides finite and causal expressions for the 
forces associated with motions of two mirrors Jl7| . 

When condition (1) is obeyed, one can ignore the recoil 
effect in the derivation of the motional force and obtain 
the susceptibility as 



XffM = im TUj 3 T[uj} 



(13a) 



where T[lo] is a cutoff function, whose expression is anal- 
ysed in detail in reference |J. A simple example, where 
the aforementioned conditions are fulfilled, is provided 
by the following lorentzian model for the reflectivity (r) 
and transmittivity (s) functions 



-1 



1 



s[u] = 1 + r[oj] 



XffM =6m rtf 

ILO 

"a " 



X 



U3 



, IUJ \ IUJ 

2Q3-1 1 - nj ln l 1_ a 



In the present paper, we will use only some general prop- 
erties obeyed by the cutoff function and not its explicit 
expression. 

The function T[u] cuts off the high frequencies and is 
regular at the low frequency limit 
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To - r[0] - rl 

ro is the reflection amplitude at w — 0; To = 1 for a 
perfectly reflecting mirror at zero frequency. It can be 
written as the sum of a dissipative and a dispersive parts 

T[u]=T R [u\+iT I {u] T R [w] and T j [w] real 

& F [u] = m TLu 3 T R [uj} (13b) 

The dissipative part T R obeys a positivity property || 

T R [u] > (14) 

The high frequency behaviour of T R [u>] will be important 
in the following. Various conditions will be expressed as 
inequalities for a positive coefficient A 



A 



d\nT R 



lim 



(15) 



Using the causal properties of the susceptibility, one can 
write a dispersion relation with subtractions since the 
susceptibility does not vanish at high frequencies. As the 
three quantities XffIP]> Xff'M an d Xff"[0] vanish, this 
relation has a simple form in terms of the cutoff function 



rM = 



<W T r [cj'} 



lit OJ 



it 



It follows that this function behaves as — at high frc- 
quencies 



«wr[cj] — > ojc for to — » oo 



(16a) 
(16b) 



Finiteness of the integral (16b) has been assumed, which 
requires A > 1. The parameter loq has the dimension of 
a frequency and may be considered as the definition of 
the reflection cutoff if To = 1; otherwise, it is the product 
of T and of the cutoff. 

In order to deal with stationary fluctuations, we will 
require that the mirror coupled to vacuum behaves as 
a stable system. Runaway solutions occur for a perfect 
mirror u>c — * oo, but not for a mirror such that 



LOcT <C 1 



(17) 



Indeed, the coupled system obeys passivity |18| in this 
case. This property is connected to the fact that vacuum 
cannot deliver energy for sustaining runaway solutions, so 
that stability follows from passivity [^) . Then, a partially 
transmitting mirror obeying equation (1) is stable. 



FLUCTUATIONS OF COUPLED VARIABLES 

We now come to the evaluation of coupled position and 
force fluctuations. 



We describe the coupling by two linear response equa- 
tions 



g[w]=9 ia H + ^[w]FH 



(18a) 
(18b) 



Each coupled variable, F or q, is the sum of the cor- 
responding input quantity, F ln or q m , and of a linear 
response to the other variable, q or F, characterized by 
the linear susceptibilities x'ff anc ^ X l q q which have been 
previously discussed. 

The linear equations (18) are easily solved, leading to 
expressions of the coupled variables in terms of the un- 
coupled ones and of linear susceptibilities which are de- 
noted Xqq, Xff, XqF = XFq for reasons which will soon 
become apparent 



q[u}=x qq [u}F in [w} + XF g [u>]q in [uj} 
F[u] - xfM^M + XffM<TM 



XqqM 
XFq[u] = 

XffH = 



X'qq M 



XffM 



1 



XfVMx^M 
i 



(19a) 
(19b) 

(19c) 
(19d) 
(19e) 



One then deduces the commutators for the coupled 
operators (£ qF = £ Fq ; ^f = & g = °) 



£z-?M = Xqq[u]CFF[u]Xqq[-u] 

+XFq[u]C q n q [u}XFq[-u} 

CFgH = Xqq[u]€FF[u]XFq[-u] 

+XF q [u]C q tt q [u]XFF[-v] 
£ FF [w] = XFq[v]gFF[v]XFq[-v] 

+Xff[w]^™[w]xff[-w] 



(20a) 
(20b) 
(20c) 



Using the fluctuation-dissipation relations (5,12) for the 
input variables, one derives, after straighforward trans- 
formations, similar relations for the coupled ones (A and 
B stand for q or F) 

2i£ A B M = Xab M - Xab [-w] (21a) 
One also obtains, after similar transformations 



C AB [uj] = 2h9(uj)UB[u} 



(21b) 
(21c) 



These expressions constitute a consistent description of 
the coupled fluctuations. The input fluctuations gr ln and 
F ln have been treated on an equal foot in the derivation 
and the resulting expressions are symmetrical in the two 
quantities q and F. The susceptibility functions Xqqi Xff 
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and \Fq play the same role as x q n q an d Xff anc ^ are 
associated to coupled variables rather than to input ones. 

We show in Appendix A that these expressions can be 
rewritten in a simpler, but less symmetric, form where 
the input position fluctuations q ln no longer appear. The 
basic reason is that the two functions XFq and xff vanish 
at the suspension eigenfrequencies ±Wo where the input 
position fluctuations g ln are concentrated. Then, equa- 
tions (19) become 



F[uj] = X Fq[uj]F in [uj} 



XqqM 
XFqM 



m (a>§ ~u) 2 ) -x™fM 



to (wq - u 2 ) Xqq[u] 
Also, the commutators (20) are simplified to 



— i,, r,.,ii 2 tfin r,.,i 



-w 



] = Xqq[u]?FF[v]XFq[-i. 
= m (cjg - lo 2 ) £ qq [u\ 

£ff[w\ = XFq[u]£F F [u\XFq[- 

= ml(LU 2 -Lj 2 ) £ qq [u] 



-UJ 



(19'a) 
(19'b) 

(19'c) 
(19'd) 



(20'a) 
(20'b) 

(20'c) 



In other words, the position fluctuations of the coupled 
mirror, as well as the coupled force fluctuations, are en- 
tirely determined by input fluctuations -F ln of vacuum ra- 
diation pressure. Although q nl and F ln have been treated 
on an equal foot, a privileged role is attributed to F ln in 
the end. As usual in the theory of irreversible phenom- 
ena, the reason for this asymmetry is that F ln corre- 
sponds to a dense spectrum whereas q m corresponds to 
a discrete one. As a result, the linear response equations 
(18) can be replaced by a simple Langevin equation (see 
Appendix A). 

This suggests that the present approach presents some 
analogy with the stochastic interpretations of quan- 
tum mechanics, like stochastic electrodynamics JL9| or 
stochastic mechanics 20 . It indeed suggests that quan- 



tum fluctuations of mirror's position can effectively be 
considered as a consequence of its coupling with zero 
point fields. As a result, these fluctuations scale as ^- 
with h introduced through vacuum fluctuations and too 
introduced as the effective mass in the Newton equation 
of motion. It has however to be stressed that the non- 
commutative character of the vacuum fluctuations plays 
a central role in our derivation, whereas the zero point 
fields are often treated as classical variables in stochastic 
interpretations. 

Moreover, it has to be noted that the reverse situation 
could occur where the field fluctuations would be domi- 
nated by position ones. Consider for example the original 
Casimir problem Km with two perfect mirrors. In this 



case, the intracavity field possesses a discrete spectrum. 
Then, one concludes that the outside field fluctuations 
determine the position fluctuations of the mirrors, turn- 
ing their discrete spectrum into a continuous one, which 
themselves determine the intracavity field fluctuations. 



COMMUTATION RELATIONS FOR MIRROR'S 
VARIABLES AND EFFECTIVE MASS 

From fluctuation-dissipation relations (21), the posi- 
tion commutator £ gg of the coupled mirror is directly con- 
nected to the mechanical admittance Y and impedance 
Z (compare with the similar relation 4 for the uncoupled 
mirror; Yr and Zr are the real parts, i.e. the dissipative 
parts, of Y and Z; a discussion of position fluctuations 
in terms of Yr and Zr is given in Appendix A) 



iuiXqq[u] =Y[w 

u£,qqM = Y R [0j] = 



1 



Zr[u>] 



(22a) 
(22b) 



In the following, we rewrite the susceptibility as (see 
eqs 13,19'c) 



1 



= toqCiJo — uj 2 m[uj] 



XqqM 

m[uj\ = mo (1 + io7Tr[(jj]) 



(23a) 
(23b) 



The coefficient too appears as the quasistatic mass, that 
is the effective mass at the low frequency limit. Using 
the relations (16) derived from the causal properties of 
the susceptibility, one deduces a different value for the 
high-frequency mass TOqo 



(X = moLOcT 



to - m 
6ttc 2 



(23c) 
(23d) 



fi appears as an 'induced mass', always positive, and is 
similar to the 'electromagnetic mass' of classical electron 
theory B (using eqs 16b and 13b) 



FF L 



(23e) 



Here, the induced mass fi is much smaller than the qua- 
sistatic mass Too for mirrors obeying equation (1), and 
low- and high-frequency values of the mass only slightly 
differ. 

Using equations (6,22), one deduces that the canoni- 
cal commutator between velocity and position is directly 
connected to the admittance function, as it is modified 
by coupling to vacuum (the uncoupled expression was the 
same with Y m in place of Y; see eq. 4) 

CM = -iu>iqq[L0] =-i(Y[u]+ Y[-u]) (24a) 
<[«(*), «(0)]> = -ih (Y(t) + Y(-t)) (24b) 
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The time dependence of the function Y{t) is analysed in 
Appendix B. Expanding in the small parameters loqt and 
ujct, assuming that loo is much smaller than u>c, we get 
the approximate expression 

r(t) = r in (t)cx P ^ + Ay(t) 

The first term corresponds to the uncoupled admittance, 
which is damped on very long times with a damping con- 
stant 

7 « WqtT 

The second part AY(t) is a 'bump' centered at t = 0, 
with a small height and a small width of the order 

° mo 

of , that is the inverse of the cutoff frequency (see eqs 
16).° 

The various frequencies involved in the admittance 
function scale as follows 



7 < w « wc < 

7 
w 



1 



(25a) 
(25b) 



Then, the time ^ is much longer than any other charac- 
teristic time 



1 f 1 

t <C — <C — < - 



(25c) 



One deduces (see Appendix B) that the commutator 
has the following values as a function of time 



<[*>(*), «(0)]> = - — 



ih l 

for t < — 

moo wc 

ih r 1 1 
for — < i < — 

m luc w 

ih , ^ ( j\t\ 

cos (u> i) exp — — 

m V 2 

for — < t 

UJ C 



(26a) 
(26b) 

(26c) 



The commutator at times shorter than the oscillation 
period, in particular, the commutator for an unbound 
mirror ujq — 0, has a time dependence directly connected 
to the difference between the two masses which appear in 
equations (23). At short times t <C and longer times 
<C t respectively, the commutator is related to the 
high frequency mass moo (mass for loq -C w) and to the 
quasistatic mass mo (mass for to <C wc)- 

This discussion suggests that moo has to be considered 
as the bare mass while m contains a part due to coupling 
with vacuum. In this spirit, we may rewrite equations 
(23) as 



m[Lo] = moo + MM 
n[u>] = itiqt (loc + iu>T[u>]) 



(23'a) 
(23'b) 
(23'c) 



where the frequency dependent function fj,[u>] appears as 
a contribution of the vacuum energy bound to the mir- 
ror, and swept along by its motion. This phenomenon 
is effective at low frequencies, but not at high frequen- 
cies since the mirror is transparent at this limit and its 
motion is decoupled from vacuum fields. 



POSITION AND VELOCITY VARIANCES 

We now come to the evaluation of the position and 
velocity variances which can be deduced as integrals of 
the symmetrical correlation function (see eqs 9). 

Before evaluating these integrals, we derive a Schwartz 
inequality for the product of the two variances 



(1 + a\io\) a qq > for any a real 

— ( x 2 

a qq Lj 2 a qq > ( \Lj\cr qq 

Using equations (21c) and (24), one gets 



\to\a qq = Lu£ qq = i£ vq (Q) 



leading to a Heisenberg inequality, which sets a lower 
bound on the product of variances, related to the equal- 
time canonical commutator 



AqAv > 



Ti 



2m c 



> 



h 

2mo 



(27) 



In order to explicitly evaluate the variances, we write 
the noise spectra as (see eqs 20') 



1 



m ( w 2 



cr qq [uj} =e(w)Z qq [w] = \£ qq [<J\ 



(28a) 
(28b) 



As t is a very small time, the noise spectrum a qq is 
approximately a sum of two narrow peaks around the 
two suspension eigenfrequencies. The width of the peaks 
is the damping constant 7 encountered previously in 
the time dependence of the commutator. Besides these 
peaks, the noise spectrum also contains a small back- 
ground (sec Appendix B). If one considers the limit where 
the mirror is decoupled from vacuum, one gets a vanish- 
ing width 7 for the peaks and a vanishing height for the 
background (compare with the discussion of the time de- 
pendence of Y(t)). Then, a simple expression is obtained, 
which corresponds to the uncoupled oscillator (see eqs 3), 
so that the usual variances (11) are recovered. However, 
this limit must be considered more carefully, as shown by 
an examination of the convergence of the integrals (9). 

Due to the high frequency behaviour of Tr[u>] (see eq. 
15), a finite value is obtained for Aq 2 as soon as A > 
(A > 1 was already needed for getting a finite value 
for toe)- It can be checked that the value for Aq 2 is 
close to the expression (11) associated with an uncoupled 
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oscillator (r — > or To — > 0) . A more stringent condition 
A > 2 is required in order to get a finite value for Av 2 ; for 
v4 = 2, a logarithmic divergence is obtained. This implies 
that the velocity variance may be completely changed 
by the coupling, even for a vanishing value of r or Tq. 
The fact that we may obtain an infinite variance for the 
instantaneous velocity is not surprising, a similar result 
being obtained in classical theory of Brownian motion. 

In the particular case of an unbound mirror (ujq = 0) 
coupled to vacuum, one gets 

£ [ w ] = — ^ 5 (29a) 

<t ot M = |£«MI (29b) 

The low frequency divergence of the spectra (29) leads 
to an infinite position variance Aq 2 , typical of the un- 
bound diffusion process experienced by a free particle 
submitted to a random force. For a mirror in vacuum, it 
could be expected that the stationary state is, as in stan- 
dard quantum mechanics, a velocity eigenstate (Aw 2 = 0; 
Aq 2 = oo). However, the integral Av 2 diverges if A < 2, 
as in the general case of a suspended mirror's motion. 



In the particular case of an unbound mirror, may 



be written (see Appendix B) 



Ay K (t) 



The contribution to the integral a' qq (t) of the first term, 
which is time independent, vanishes. Since AYr differs 
from zero only at short times, of the order of — , the 
asymptotic behaviour at long times of a' qq (t) can be eval- 
uated as 



,(*) 



A1«[0] 



Tft 



tT 
irmot 



for uic<t ^> 1 



It follows that the quantum diffusion of an unbound mir- 
ror coupled to vacuum is characterized by a logarithmic 
behaviour at long times E4] 



A(t) 



UtTq t 
In — 

7rm t 



for uoct ^> 1 



to is an integration constant, of the order of 



This 



diffusion corresponds to very small displacements, the 
length scale being in fact the Compton wavelength 
associated with the mirror 



QUANTUM DIFFUSION OF A MIRROR 
COUPLED TO VACUUM 

The quantum diffusion of a mirror coupled to vacuum 
may be characterized quantitatively by the function 

A(t) = \ ((q(t) q(Q)) 2 ) = h (a qq (Q) - a qq {t)) 



2^ 



ha qq [ui] (1 — cos (u>t)) 



This integral is well defined for any value of t, even for 
an unbound mirror (ujq = 0), and it vanishes when t goes 
to zero. This is also the case for the integral associated 
with A'(t) (we suppose A > 1), but not for A"(t), which 
is not defined if A < 2 (A"(0) = Av 2 ). 

Interesting results about quantum diffusion are ob- 
tained by using the analytic properties of the correlation 
function C qq (t). Since the noise spectrum C gg [w] con- 
tains only positive frequency components (see eq. 21b), 
C a q (t) is analytic and regular in the half plane Imt < 
p2| , f23|| . Then, <J qq {t) may be deduced from £, qq (t) through 
a dispersion relation (V stands for a principal value) 



a qq {t) = dt'£ qq (t-t')V 



1 

iTTt 7 



(30) 



One obtains a simple integral relation between the time 
dependent commutator and the position variance by set- 
ting t — in the preceding integral. One also obtains by 
deriving the expression of a qq (t) 



dt'e Q M-t')v 



i 

iirt' 



A(t) 



6tt 2 



moc 



In (tuct) 



(31) 



It has to be noted that thermal effects [[15] change the 
behaviour of A(t) at long times, even at a low tempera- 
ture. 



FLUCTUATIONS OF COUPLED FORCE 

It is interesting to discuss the autocorrelation of the 
coupled force F, that is the force computed in presence 
of the mirror's radiation reaction, and its correlation with 
the coupled position. 

The autocorrelation of the coupled force F is given 
by equation (20'c); £ff is identical to £p F , except in 
the narrow resonance peaks where the fluctuations of the 
coupled force are smaller than those of the input force 
(compare with eqs 28) 



Zff[w]=&f[u] 



(- 2 -- 2 r 



LOP. 



(32) 



The correlation of the coupled force with the coupled 
position is described by equation (20'b). This expression 
implies that coupled position and coupled force evaluated 
at equal times commute; is an odd function of u 

as C 9g M, so that 



Wo) = o 



(33) 



From the same argument of parity, it follows that coupled 
velocity and coupled force do not necessarily commute 
when evaluated at equal time. 
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One also deduces the symmetrical correlation function 
between q and F 



The estimator q[uj] for a frequency component q[u] of 
position can be written 



a Fq [uj] = e(w)£ Fq [w] = m (u 
(g(O)F(O)) = (F(O)g(O)) = fia^ 
= mo (woAg 2 - Av 2 ) 



'qqi 



(34) 



which is reminiscent of the virial theorem |2q| . 

It appears from this discussion that the position of the 
mirror and the vacuum pressure exerted upon it are in- 
timately correlated: at each frequency, they are directly 
proportional to each other. Therefore the quantum fluc- 
tuations of the mirror cannot be treated independently 
of those of vacuum radiation pressure. 



2k 

The added noise gjvM is a sum of two terms associated 
respectively with phase and intensity fluctuations (com- 
pare with eqs 11 of ref. ||) 

5$[u)] 

QnM = -rpr H 2UkaXqq[u]5I[(d] 

ZKq 

When phase and intensity fluctuations are considered 
as uncorrelated, the corresponding noise spectrum is ob- 
tained as a sum of two contributions lEI 



ULTIMATE QUANTUM LIMITS IN A POSITION 
MEASUREMENT 

We show in this section that the noise spectrum a qq 
can be considered as the ultimate quantum limit for the 
sensitivity in a position measurement, when mirrors cou- 
pled to vacuum are used. More precisely, we discuss how 
this limit may be reached, at least in principle; we will 
not analyse in detail the effective realization of this sen- 
sitivity. 

We consider a simple model of position measurement 
where a laser is sent onto the mirror and where the phase 
$ of the reflected beam is monitored || . This phase $ is 
related to the mirror's position q through (feo is the laser 
wavenumber) 

= 2k {q(t) + 5q RP (t)) + 5$(t) 

Here, q(t) represents the position discussed in foregoing 
sections, with its quantum fluctuations; <5<£> represents the 
quantum fluctuations of the phase of the incident field, 
and Sq R p stands for the random motion of the mirror 
due to the radiation pressure 5F RR exerted upon it |Q, 
proportional to the fluctuations SI of laser intensity 

SF RP {t) = 2hk SI(t) 

In the following, we will consider that the fluctuations 
of q, which can be assigned to vacuum radiation pres- 
sure (the vacuum force at frequency ui is associated with 
vacuum field fluctuations at frequencies comprised be- 
tween and uj), are not correlated with the fluctuations of 
8q R p and (field fluctuations at frequencies close to the 
laser frequency). We will also assume that the field fluc- 
tuations around the laser frequency can be transformed 
(squeezed) at our convenience, while the lower frequen- 
cies vacuum fluctuations are not modified. 

In a linear analysis, the response of the mirror to radi- 
ation pressure is described by the susceptibility function 

SqRp[(d] = Xqq[u]8F RP [uj] 



tTtfcj) [id] 
(2fco) 2 



(2hk ) 2 \x qq [td}\ 2 <Jii[u;} 



The minimum value reached by this expression when the 
various parameters are varied, with variations of o$$ and 
an constrained by a Heisenberg inequality, is the stan- 
dard quantum limit 



It is known that the sensitivity can be pushed beyond 
the standard quantum limit Q . Indeed, phase and inten- 
sity fluctuations are linearly superimposed in the fluctua- 
tions of the monitored signal, which therefore depend also 
upon the correlation 07$. Consequently, one can reduce 
noise by squeezing a well chosen quadrature component 
of the fields. A lower bound ct u Q l is obtained by opti- 
mising the parameters characterizing the squeezed fields, 
the variations of 0$$, an and 07$ being constrained by a 
Heisenberg inequality. It is determined by the dissipative 
part of the mirror's susceptibility || . With notations of 
the present paper, this 'ultimate quantum limit' can be 
written 



We have used relation (28b) between the noise spectrum 
and the dissipative part of the susceptibility. 

It is instructive to measure the fluctuations of the esti- 
mated position gas a noise energy N per unit bandwidth 
(fluctuations of q and q^ are uncorrelated) 



JVM = 



to 



(C, 



qql 



a 



qnqN I 



A) 



The standard quantum limit corresponds to 

JV SQL M = h9(td)m a x 

H+u 2 ) (KmMI + IxtoMI) 



ne {id) 



for ujq <C (d <C id S 



(35a) 
(35b) 



This limit is entirely determined by added fluctuations 
(second term in 35a) and proper quantum fluctuations of 
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position (first term in 35a) do not contribute (these fluc- 
tuations are concentrated in a very narrow peak at the 
suspension eigenfrequency) . This confirms that proper 
fluctuations can be forgotten when analyzing a position 
measurement with a sensitivity level around standard 
quantum limit. 

In contrast, the ultimate quantum limit reaches the 
level set by the dissipative character of coupling with 
vacuum; proper fluctuations have the same contribution 
as added fluctuations 

N vql [lu] = h9(Lo)m a (u% + lu 2 ) 2\£ qq [u}\ (36a) 
« 2M(u)u)tT for uj q <C ui <C uj s (36b) 

This limit (36) is far beyond the standard quantum limit 
(35). It has to be stressed however that the input state 
of the field used in the measurement must be carefully 
controlled, in the frequency band around laser frequency 
which is involved in the measurement. Reaching effec- 
tively the ultimate limit in a real experiment may be a 
rather difficult task. 

This is made clear by expressing this limit as a position 
variance 

Aq 2 a r^Hii ( — Y 
u> 3ir \moc J 

for wq Cw - B,uj + B <^ lus 

where is the Compton wavelength of the mirror 
and 2B the measurement bandwidth (signal measured 
through a filter characterized by a function G[u>] having 
a maximum value of 1 at expected signal frequency) 

2B = G 

The effective noise Aq 2 in the measurement may be much 
smaller than the position variance: the integral of a 
positive quantity over a limited bandwidth is obviously 
smaller than the integral of the same quantity over all 
frequencies (see eqs 9). This occurs, if the input state of 
the laser is optimized, when the suspension eigenfrequen- 
cies are outside the detection bandwidth; in the opposite 
case, the effective noise would be of the order of the po- 
sition variance. Incidentally, this discussion means that 
quantum noise may be reduced by a narrow band detec- 
tion, exactly in the same manner as thermal noise pq ], 

DISCUSSION 

Any mirror, more generally any scatterer, is coupled to 
vacuum radiation pressure, so that its position acquires 
quantum fluctuations from vacuum. This mechanism is 
quite analogous to the thermalization of the variables of a 
system coupled to a thermal bath, but it involves vacuum 
(i.e. zero temperature) fluctuations. 

Quite generally, coupling a classical quantity to quan- 
tum ones leads to inconsistencies j27|. In contrast, 



generic quantum properties are preserved in a consistent 
treatment, as a consequence of fluctuation-dissipation re- 
lations. Then, position fluctuations of the coupled mirror 
are connected to the dissipative part of its mechanical ad- 
mittance. 

It appears that fluctuations of the coupled variables 
are determined by input fluctuations corresponding to a 
continuous spectrum (i.e. the 'dissipative' ones), which 
are here the force fluctuations. It can therefore be as- 
serted that quantum properties of position are generated 
from vacuum fluctuations: the same results would be ob- 
tained if the position were initially treated as a classical 
number. Nevertheless, the proper position fluctuations 
are recovered in the end, in the limit of decoupling. Cou- 
pling with vacuum is very weak, making the value of r 
extremely small, and this justifies that usual quantum 
mechanics can provide us with a good description of po- 
sition fluctuations. 

The canonical commutation relation between position 
and velocity is slightly modified, this modification being 
only effective at small times of the order of the reflection 
delay. This time dependence is connected to the differ- 
ence between the bare mass (at high frequencies where 
field energy cannot follow the mirror's motion) and the 
quasistatic mass (at low frequencies where a finite part of 
the vacuum energy is swept along the mirror's motion). 

The noise spectrum, which describes the quantum fluc- 
tuations of position for a mirror coupled to vacuum, con- 
sists in two narrow resonance peaks close to the sus- 
pension eigenfrequencies and a small broad background. 
When the peaks are approximated as Dirac distributions 
and the background disregarded, integrals over frequency 
of the noise spectrum reproduce the dispersions Aq and 
Av expected for a quantum harmonic oscillator in its 
ground state. 

However, differences exist which may have some impor- 
tance. Although finite in the uncoupled case, the veloc- 
ity dispersion may be infinite for the coupled oscillator, 
a property reminiscent of classical theory of Brownian 
motion. Also, the resonance peaks have a small width, 
typical of dissipative coupling with vacuum, so that the 
correlation functions are damped on very long times and 
differ from the usual dissipationless expressions. For an 
unbound mirror coupled to vacuum, a logarithmic diffu- 
sion is obtained (no diffusion for an unbound uncoupled 
mirror), the length scale of which is given by the Comp- 
ton wavelength associated with the mirror. 

Finally, the proper quantum fluctuations of mirrors 
used in a length measurement such as interferometric de- 
tection of gravitational waves are accounted for as soon 
as coupling with vacuum radiation pressure is treated. 
The noise spectrum associated with position sets the ul- 
timate bound on sensitivity when the measurement is 
optimized. Then, the effective noise in the measurement 
can be much smaller than the position variance computed 
from quantum mechanics, as soon as the detection band- 
width does not contain the suspension eigenfrequencies. 
Quantum noise, as other sources of noise, may be reduced 
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by monitoring a signal away from resonance frequencies. 
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APPENDIX A: ELIMINATION OF THE INPUT 
POSITION FLUCTUATIONS 

The linear response equations (18) can be replaced by 
simpler, but less symmetric, ones where the input posi- 
tion fluctuations q m no longer appear. 

The demonstration is based upon the following rela- 
tions (deduced from eqs 19) 



XFql 



Xqql 



= mo (wq -UJ 2 ) Xqq[w] 



XffM = Xf 9 MxffM = m o ( w o - w 2 ) XwHXfVH 

It follows that the functions XFq and xff vanish at the 
suspension eigenfrequencies 

XFq[±u ] = xff[±uj ] = 

Since the input position fluctuations q ln are concentrated 
at these frequencies (see eqs 3), they cannot feed the fluc- 
tuations of the coupled variables which are therefore en- 
tirely determined by the input force fluctuations. Equa- 
tions (19) become 

q[u>}=xMF in [uj] 
F[w] = to (w 2 - uj 2 ) q[uj] 

Xqq[u] 



too (uj 2 - uj 2 ) - XffM 



Although q ln and F ln have been treated on an equal foot, 
F ln plays a dominant role in the end. This is due to the 
property 



XqqM 



= 



which is satisfied for the discrete spectrum of equations 
(3). The same property is not fulfilled by the dense spec- 
trum associated with the force (see eq. 13). It results 
that expressions (20) of the commutators can be rewrit- 
ten as (20')- 

A mirror coupled to vacuum then obeys a simple equa- 
tion of motion where, as in Langevin theory of Brownian 
motion £3], the force F is a sum of a Langevin force F m , 
the input force, and of a long term cumulative effect, the 
motional force, 

m 0J 2 q(t) + m q"(t) = F(t) 
to (uj 2 - uj 2 ) q[u>] = F[u>] 
F[u] = F in [u] + XffIvMu] 



These equations can be obtained directly by setting q ln = 
in (18). 

These equations can be solved in terms of a mechanical 
impedance Z or of a mechanical admittance Y, which 
connect the applied force and the effective velocity 

F in [uj] = Z[uj]v[uj] 

Z[uj] = — - — - — iujrriQ (1 + ic<jTrL>]) 

— iuj 

v[uj] = Y[uj]F in [uj] 

The dissipative parts of the impedance and admittance 
functions, defined as their real parts Zr and Yr in the 
frequency domain, are positive at all frequencies (see eq 
14) 



Zr[uj] = m Q uj 2 TT R [uj} > Y r [uj] = 



Zr[uj] 



> 



In fact, the impedance and admittance functions are pas- 
sive functions, which ensures stability of the coupled sys- 
tern §. 

We can then give a simpler derivation, in the spirit 
of Nyquist's derivation O-O], of the fluctuation- 
dissipation relations for coupled position. We first 
rewrite the relation between the commutator and 
the dissipative part of the motional susceptibility Xff 
(see eqs 12) 

& F [w] = ujZ r [uj] 

Hence, the non commutative character of velocity fluc- 
tuations is related to the dissipative part Yr of the me- 
chanical admittance 

= I^MI 2 £ffM = uY R [uj\ 

The result for the position's commutator is the same as 
previously obtained (see eqs 21) 

£™M Yr[uj] 

ZqqM - - — 



APPENDIX B: TIME DEPENDENCE OF THE 
COMMUTATOR 

In order to analyze the time dependence of the coupled 
commutator ^g[cj], we first consider the case of an un- 
bound mirror coupled to vacuum, whose admittance can 
be written (luq = 0) 



Y[u] = 
AY[uj] = 



1 



1 



(e — iuj) m[uj] too (e — iuj) 
toq (1 + iwrTfw]) 



AY[uj] 



Y(t) = M + AY(t) 
mo 
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In this decomposition, the first part is the uncoupled ad- 
mittance while the second part represents a correction. 

One deduces Y(0 + ) from the high frequency behaviour 
of Y[u] 

Y(0+) = lim (-iujYU) = — 

Starting from the value for t = + , Y(t) decreases 
to the slightly different value for t — > oo. In other 
words, Y(t) is the sum of a constant value — and of a 

' \ > mo 

bump AY(t), having a height 

111 U>C'T u>c T 

m 



One then deduces by straightforward calculations 



AF(0+) = 



m 



moo mo too 1 - w c t 
The integral of the bump is given by 

AF[0] = ^ 

m 

so that the width of the bump is of the order of ^ SL . 
One deduces from equations (24) 

if) 

<[«(*), «(0)]) = ih (AY(t) + AY(—t)) 

m 

The first term is a canonical commutator (see eq. 10) 
while the second is a bump-shaped correction, having a 
width of the order of the mean reflexion delay — . At 

•> LUC 

short and long times, one finds equations (26a,b) 



(Ht),q(0)]) 



ih r 1 
for t <C 

moo w c 

in r 1 

for — < i 

mo wc 



For an harmonically suspended mirror, the admittance 
function Y can be written as a sum of three components 
rather than two 

Y[lu] = . P+ . + , P " , + AY[uj] 

Y(t) = 6(t) ( P+ exp(-iuj + t) + p_ exp(-iw_t)) + A7(t) 

The two complex numbers w± are the poles of the admit- 
tance function Y[u], which are close to the suspension 
eigenfrequencies ±ojq 



P± = 



mo^Q = aj±m[w±] 



2m[w±] + w±m'[u;±] 



We can estimate these quantities in an expansion with 
respect to the small parameters uqt and u>ct. Assuming 
that uoq is much smaller than wc, we get at lowest order 

7 « cjgTr 



o 2 T 



AY[0] 



Tr w c t 



too 



AF(0+) 



m 



so that the shape of the function Ay is similar to that 
previously obtained for an unbound mirror. 

Consequently, at short and long times, one finds 
(26a,c) 

([«(*), g(0)]> =-■£- fort«-l 

TOqo 

TO ,_ . / j\t\ 

= cos (uj t) exp — 

m V 2 

for — <C t 

0J C 



APPENDIX C: CASE OF AN ANHARMONIC 
SUSPENSION 

For the sake of simplicity, the simple case of an 
harmonically suspended mirror has been considered 
throughout the paper. We show here how the results can 
be generalized to the case of an anharmonic suspension. 

Denoting uj a the various excitation frequencies from 
the ground state to the upper states, and q a the corre- 
sponding matrix elements of the position operator, one 
obtains 

C£(*) = X>a| 2 exp(-iu; Q t) 



Hence 

4" H = ^ E \Q»\ 2 ** (*(<■> <■>«) - + "a)) 

a 

c£M - E kal 22 ^ - w «) = 2 ^)C M 

<M = ^^|<Za| 2 27r(%-w a )+<5(u; + u; a )) 

a 

One also defines the susceptibility 



2w n 



ft^'" 1 "' w?, + (e — iui) 2 

x£M-x£[-"]=2i^[a;] 

The frequency dependence of the susceptibility is 
less simple than in the harmonic case. However, the 
fluctuation-dissipation relations remain valid, for coupled 
as well as uncoupled variables. Also, the following prop- 
erty remains true for a discrete spectrum 



P± 



m (2 + 3iu>±TT ) 
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1 

: r~ r = 

X&[±Wa] 

QH =0 

which allows one to express the coupled fluctuations in 
terms of the uncoupled force fluctuations only. 
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